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The epoch of the end of reionization and the Thomson optical depth to the cosmic microwave background 
depend on the power spectrum amplitude on small scales and on the ionizing photon emissivity per unit mass in 
collapsed halos. In this paper is investigated the role of the radiative feedback effects for the temporal evolution 
of the ionizing emissivity. 

It is shown that the observational constrains on hydrogen photo-ionization rate based on Ly-a, Ly-/3 and Ly-7 
Gunn-Peterson troughs and an electron optical depth consistent with the latest CMB measurements requires an 
emissivity of ~10 ionizing photons per baryon and Hubble time at 2 = 6. 

Through E-mode CMB polarization power spectrum measurements, is expected that Planck experiment will 
have the sensitivity to distinguish between different histories of stellar evolution. 



1. Introduction 

The detailed study of the intergalactic medium 
(IGM) is fundamentally important for under- 
standing the large-scale structure properties and 
the galaxy formation process. 
At the epoch of reionization the collapsed objects 
began to influence the diffuse gas in the IGM 
rendered it transparent to the ultraviolet pho- 
tons. In order to virialize in the potential wells 
of the dark matter halos, the gas in the IGM 
must have a mass greater than the Jeans mass 
Mj ~ 1O 5 M0 at a redshift z ~ 30, corresponding 
to a virial temperature of T V i r ~ 10 4 K. Photoion- 
ization by the high-redshift UV radiation back- 
ground (UVB) heats the low density gas in the 
IGM before it falls into the dark matter wells, 
strongly reducing the fraction of neutral hydro- 
gen and helium that dominate the cooling of the 
primordial gas at temperatures of T V i r . 

Reionization is an inhomogeneous process that 
proceeds in a patchy way. The radiation out- 
put associated with the collapsed halos gradually 
builds up a cosmic UVB. At early times, most of 
the gas in the IGM is still neutral and the cosmo- 
logical Hji regions around the individual sources 
do not overlap. At this early stage the gas in the 
IGM is opaque to the ionization photons, caus- 



ing fluctuations of both ionization fraction and 
UVB intensity from region to region. At the 
reionization redshift the Hu regions surrounding 
the individual sources in the IGM overlap. Dur- 
ing reionization, the number of ionizing photons 
reaching IGM had to be sufficient to ionize every 
atom in the universe and to balance the recombi- 
nation. 

From the observational point of view, the study 
of the reionization process itself as well as the 
properties of the sources driving it is challenged 
by a variety of observational probes. A powerful 
observational probe comes from the Lya absorb- 
tions spectra of the high redshift quasars (Becker 
et al., 2001, Fan et al., 2002, White et al, 2003, 
Fan et al., 2004) showing that all known quasars 
with z > 6 have a complete Gunn-Peterson (GP) 
trough and a rapid evolving hydrogen neutral 
fraction compatible with the final stage of reion- 
ization. An other powerful observational probe 
is represented by the high value of the electron 
scattering optical depth r inferred from the Cos- 
mic Microwave Background (CMB) anisotropy 
measured by the WMAP experiment (Kogut et 
al., 2003; Spergel et al., 2003; Verde et al., 2003) 
which requires reionization to begin at z > 14. 

In this paper we study the temporal evolu- 
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tion of the ionizing emissivity in the presence 
of the radiative transfer effects by using N-body 
cosmological hydrodinamical simulations at sub- 
galactic scales. Detailed discussion of the simula- 
tions can be found in Popa, Burigana, Mandolesi, 
2005 (hereafter PBM). Throughout it is assumed 
a background cosmology consistent with the most 
recent cosmological measurements (Spergel et al., 
2003) with energy density of fl m = 0.27 in mat- 
ter, Ob = 0.044 in baryons, Oa = 0.73 in cos- 
mological constant, a Hubble constant of H =72 
km s _1 Mpc _1 , an rms amplitude of o% = 0.84 for 
mass density fluctuations in a sphere of radius 
8h _1 Mpc, adiabatic initial conditions and a pri- 
mordial power spectrum with a power-law scalar 
spectral index n s = 1. 

2. Radiative transfer 

The cumulative UV background flux, J{v , z D ), 
observed at the frequency v Q and redshift z Q , (in 
units of 10~ 21 erg cm~ 2 s _1 sr _1 ) due to pho- 
tons emitted from redshifts between z a and an 
effective emission screen located at z sc > z Q , is 
the solution of the cosmological radiative transfer 
equation (Peebles, 1993; Haiman, Rces & Loeb, 
f997): 
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where: j(v z ,z) is the comoving emission coef- 
ficient (in units of 10~ 21 erg cm~ 3 s _1 sr _1 ) 
computed at emission redshift z and photon fre- 
quency v z = v Q (l + z)/(l + z ), T e ff(v ,z ,z) 
is the effective optical depth at the frequency 
v due to the absorption of the residual gas in 
the IGM between z a and z and (dt/dz)^ 1 = 
—Hq(1 + z) ^/O m (f + z) 3 + Q\ is the line element 
in the ACDM cosmology. 

Above the hydrogen ionization threshold the UV 
radiation background is processed due to the ab- 
sorption of residual gas in the IGM dominated by 
neutral hydrogen and helium. At these frequen- 
cies the effective optical depth is given by (see e.g. 
Haiman, Rees & Loeb, 1997): 



T ef f(v ,z ,z) = c J ^n(v z ,z)dz 
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Figure 1. Left panel: Evolution with redshift of 
the UVB flux at wl Ryd. The critical UVB flux 
below which (T vir < I0 2 4 K) the star formation 
is prevented (thick dotted line) is from Haiman, 
Abel & Rees (2000). Right panel: Redshift evolu- 
tion of the hydrogen photo-ionization rate for the 
UVB models presented in the right panel. The 
observational upper limits based on Lya, Ly/3 and 
Ly7 Gunn-Peterson troughs at z = 6.05 are from 
Fan et al. (2002). 



k{v z ,z) = ^2 u i{ v z)ni{z) , 

i 

where: i = (Hi , He{) and crj and n>i are the cross 
section and the number density respectively, cor- 
responding to species i. 

The hydrogen photo-ionization rate, r 7fJ ^(z), 
is related to the UVB flux, J(v, z), through: 

J(v,z)a Hl {v) 
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where ohi(v) is the hydrogen photo-ionization 
(dissociation) cross-section (Abel et al. 1997) and 
VH[ is the hydrogen ionizing threshold frequency 
(v Hl wlRyd). 

Figure 1 presents (left panel) the evolution with 
redshift of the UVB flux at « IRyd obtained 
for different UVB spectrum models (PBM) that 
constraint the hydrogen photo-ionization rate at 
z ~ 6 to T 1HI ~ 8xI0~ 14 s _1 (right panel) as in- 
dicated by the observational upper limits based 
on Lya, Ly/3 and Ly7 Gunn-Peterson troughs 
(Fan et al. 2002). In Figure 2 are presented 
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Figure 2. The evolution with redshift of the ion- 
ization fraction obtained for the UVB models pre- 
sented in Figure 1 with electron optical depths 
T es =0M - 0.1 (PBM). 



the corresponding ionization histories with opti- 
cal depths in the range r es =0.05 - 0.1 (PBM). 
The main effect of the UV radiation spectrum on 
the temporal evolution of the ionization fraction 
is given by the value of the reionization redshift 
and the redshift interval in which the reioniza- 
tion is completed. As reionization proceeds, the 
radiative feedback raises the Jeans mass in the 
ionized regions, suppressing the formation of low 
mass systems. The net effect is a decrease of the 
global ionization fraction with cosmic time over a 
limited period. 

3. Emissivity evolution 

The hydrogen photo-ionization rate is related 
to the ionizing emissivity, e# r (v, z) through: 

r-OO 

r 7ff J ( z )=/ dve Hl {v,z),\{v,z)(TH I {v), (4) 

where z) is the mean free path of the ionizing 
photons. Assuming cr# 7 [y) ~ i^ 3 , then \(v, z) can 
be written as (Miralda-Escude, 2003): 

A(i/,*) = Ao( — ) , (5) 

where Aq 1 (z), the absorbtions probability per 
unit length for a photon at v = vhi, is given 



by: 

Ao 1 W = £, ALL = cff- 1 (z)(^)" 1 I (6) 

with h{z) = H a ^n m {i + zY + n v . 

In the above equation dN^L/dz is the Lyman 
limit system abundance per unit redshift interval. 
Storrie-Lombardi et al. (1994) found dN^L/dz = 
3.3 ± 0.6 at z=4 while Madau, Haardt & Rees 
(1999) report a value 1.5 times larger (see also 
Miralda-Escude, 2003 and references therein). 

The mean ionizing emissivity integrated over 
frequency, expressed as the number of ionizing 
photons per baryon and Hubble time, can be writ- 
ten as: 

<z) = V^dNLL f°° J(y,z) (JL.Y 1 ' 6 dv 
H(z)n b c dz J VHi hp \u Hl J 

where n b = 2.07 x 10- 7 (f7 b /i 2 /0.022)(l + z) 3 Cm - 3 
is the total number density of baryons. 
Figure 3 presents the evolution with redshift of 
the mean ionizing emissivity obtained for various 
reionization histories. 

Figure 4 presents the dependence of the mean 
ionizing emissivity on the circular velocity and 
on the virial mass. For all cases, the experimen- 
tal constraint on hydrogen photo-ionization rate 
predicts ~10 ionizing photons per baryon and 
Hubble time at at z = 6. The electron optical 
depth of r es =0.05 - 0.1 requires a mean ionizing 
emissivity of earlier low-mass halos (virial mass 
of <~ 10 6 M©) of order unity. 
The emissivity of halos at z =6 (v C i rc «50 kms -1 ) 
is ^35 photons as obtained by dividing the mean 
emissivity of 10 photons per baryon and Hub- 
ble time by the baryon collapse fraction (see Fig- 
ure 5): 

F b (z) - x e f coU (z,T h ) + (1 - x e )fcou(z,T c ) , (8) 

where f C oll(z,T h ) and f co ii(z,T c ) arc the col- 
lapse fractions for the regions with the tem- 
perature higher than the virial temperatures 
T h = 2.5 x 10 4 K and T c = 10 2 4 K (Haiman, Abel 
& Rees, 2000; Barkana & Locb, 2001). 
In Figure 5 these numbers are compared with the 
number of photons that can be produced by star 
formation assuming that iV 7 =2.5xl0 4 photons 
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Figure 3. Evolution with redshift of the mean 
ionizing emissivity (number of ionizing photons 
per baryon and Hubble time) corresponding to 
the reionization histories presented in Figure 2. 



per baryon are emitted by zero-metalicity stars 
with M > 20M Q (Onken & Miralda-Escude, 
2004) and dN LL /dz = 3.6 (Storric-Lombardi et 
al. 1994). If /* is the fraction of the gas in halos 
that can be converted into stars with M > 20M Q 
and f esc is the fraction of the photons that escape 
into IGM, then the experimental constraint on 
r 7HJ ~ 8xl0- 14 s _1 requires / esc /* w 10" 19 at 
z = 6. 

The electron optical depth r es =0.05 - 0.1 pre- 
dicts Use]* = ICT 2 ' 7 - 10" 1 ' 9 for z =15 - 6. 
This values comparable with f esc f* >(10~ 2 - 4 ,10~ 2 - 8 ) 
found by Onken & Miralda-Escude (2004) for the 
same population of stars by using an analytical 
model normalized to 7 photons per baryon and 
Hubble time at z = 4 and requiring r es =0.17 
(Miralda-Escude 2003). 



4. Implications for the CMB E-mode po- 
larization 

The CMB anisotropy temperature, Ct(1), and 
E-mode polarization, Ce(1), power spectra cor- 
responding to the various ionization histories are 
presented in Figure 6. 

While Ct(1) power spectra are almost degener- 
ated, the differences in different reionization his- 
tories produce undegenerated signatures on Ce{1) 
power spectra at low multipolcs (I < 50). This 



Figure 4. Dependence on circular velocity (virial 
mass) of the mean ionizing emissivity correspond- 
ing to the reionization histories presented in Fig- 
ure 2. 
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Figure 5. Fraction of mass collapsed into halos as 
function of redshift. 
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Figure 6. Evolution with redshift and IGM tem- 
perature of the product / esc /* corresponding to 
the reionization histories presented in Figure 2. 
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Figure 7. The CMB angular power spectra, Ct(1) 
and Ce(1) for the reionization histories presented 
in Figure 2. 
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can be explained by the fact that while polariza- 
tion is projecting from the epoch of reionization 
at angular frequencies I = k(ij — r) ri ) ( here k is 
the wave number, 770 and r] r i are the conformal 
times today and at the epoch of reionization) the 
temperature is projecting from the (further) last 
scattering surface. 

It is expected that through E-mode CMB po- 
larization power spectrum measurements, the 
Planck experiment will have the sensitivity to 
distinguish between different histories of stellar 
evolution even they imply the same optical depth 
to electron scattering and degenerated Ct power 
spectra (PBM). 
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